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SUMMARY 


A linear small- deflection theory is developed for the elastic 
behavior of inflatable plates of which Airmat is an example. Included 
in the theory are the effects of a small linear taper in the depth of 
the plate. Solutions are presented for some simple problems in the 
lateral deflection and vibration of constant-depth rectangular inflat- 
able plates. 


INTRODUCTION 


For certain types of satellite and reentry vehicles it is desirable 
to have a structure which can be packaged in a compact form for launching 
and erected after injection into orbit or at the time of reentry. One 
possible way of meeting this requirement is to utilize an inflatable 
structure. Inflatable structures also have a variety of other applica- 
tions where it is desirable to have a small package for transporting and 
simple erection capability at the destination. These applications range 
from inflatable airplanes to erectable living quarters. Where platelike 
structural components are needed in an inflatable structure, for example, 
for the lifting surfaces, fins, and control surfaces of a reentry glider, 
an inflatable plate such as Airmat (developed by Goodyear Aircraft 
Corporation) appears to be an efficient and useful scheme. This type 
of plate is illustrated schematically in figure 1. It consists of two 
woven covers having airtight coatings and held some distance apart by 
the combined action of drop cords and internal pressure. The drop cords 
are closely spaced, distributed throughout the plate, and may be of 
varying length to form a plate of variable thickness. In a reentry- 
glider application, for instance, metal wire would most likely be used 
to weave the covers and form the drop cords. 

In order to make rational stress, deflection, and aeroelastic anal- 
yses of inflatable -plate structures, it is necessary to have a theory 
from which stresses and deflections in such plates can be calculated for 
various external loading conditions. In this report a linear theory for 
inflatable plates is derived, and some solutions to elementary static 
deflection and vibration problems are presented. 
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SYMBOLS 


11 ' 


a,b 


E W t 

Ept 


^12* ^21*^22*^33 orthotropic elastic constants defined in equa- 
tions ( 13 ) 

length of rectangular plate in x- and y- directions, 
respectively 

root chord length of triangular plate 
extensional stiffness in warp lirection of cover 

extensional stiffness in fill direction of cover 


e 

Gt 

h 

l 




M. 


i x , My, M X y 


unit deformation through thickness of plate 
shear stiffness of cover 
depth of plate 

semi span of triangular plate 
plate moment resultants 

external moments applied to edges of rectangular plate, 
positive when causing tension or positive shear in upper 
cover 


M Nx ,M N'\ 


N N N 

A V V xy 


N° N° N° 
X> A y j A1| xy 


vector components of moments applied to edge of plate, 

positive in the positive x- and y-directions, respectively 

number of half waves in x- and y-directions, respectively 

stress resultants associated with xyz coordinate system 
and dependent on displacements 

stress resultants associated with xyz coordinate system 
and independent of displacements 


N x N x + N x 

- *2 + N y 


N* = N° + N 
xy xy xy 


J (J\H Q\ 
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N^y stresses in covers of inflatable plate associated with 
x^’z 1 coordinate systems 

^Nx>^Ny vector components of middle plane forces applied to edge 

of plate , positive in the positive x- and y-directions, 
respectively 

p internal pressure 

Q-x^yj^z, x ~ > y - } an( i z-components, respectively, of external dis- 

tributed loading per unit middle-plane area 

S^, S transverse shear carried in covers due to taper (see eqs. (22)) 

s coordinate along edge of plate (see fig. 1(a)) 

t thickness of each cover 

u, v, w displacements in x- , y-, and z-directions, respectively 

u',v',w' displacements in x’ - , y ’ - , and z '-directions, respectively 

AV change in volume of inflatable plate 

external lateral loads applied to edges of rectangular 
plate 

V N external lateral load applied to edge of plate of arbitrary 

shape, positive in the positive z-direction 

W work of external loads and body forces 

x, y, z rectangular Cartesian coordinate system 

x',y',z T coordinate systems associated with covers 

a, 3 angles of rotation of the drop cords in the xz- and yz- 

planes, respectively, from their initial position normal 
to the xy- plane 

e',e' r , y' strains in covers of inflatable plate 

x y xy 

5 variational operator 

rectangular Cartesian coordinates locating the final posi- 
tions of points on the surfaces of the inflatable plate 
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A 

^WF 


M-fw 

n 

n c 

% 

p 

O) 


length- width ratio of rectangular plate, b/a 

Poisson’s ratio associated with a contraction in the fill 
direction caused by a tensile stress in the warp direction 

Poisson’s ratio associated with a contraction in the warp 
direction caused by a tensile stress in the fill direction 

total potential energy 

strain energy in covers of inflatable plate 
work done against internal pressure 

mass of plate per unit middle-plane area 
circular frequency of vibration 


Subscripts: 

F,W pertain to fill and warp directions, respectively 

m, n integers 

+,- pertain to upper and lower covers, respectively 

A comma followed by a subscript denotes partial differentiation with 
respect to the subscript, A dot over a symool denotes partial differen- 
tiation with respect to time. 


BASIC ASSUMPTIONS AND COORDINATE SYSTEMS 


The inflatable plate shown in figure l(a) is of arbitrary shape in 
planform and may have a small linear taper in depth in the x- and y- 
directions. The taper is assumed to be symmetric about the middle sur- 
face. The covers of the plate are assumed to be identical and are 
treated as orthotropic membranes. The closely spaced drop cords are 
assumed to be straight and inextensional, and they are conceptually 
spread continuously over the plate. The drop cords are assumed to be 
normal to the middle surface before deformation, but they are assumed 
to be hinged at the ends so that during deformation the angles between 
the drop cords and the middle surface or the covers may change. The 
internal pressure in the plate is assumed to be constant during defor- 
mation. Sidewalls are presumed to be present at the edge of the plate 
to contain the pressure, but their effects on plate behavior are not 
otherwise taken into account. 


mo h vo f 
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In this investigation the undeformed shape of the plate is assumed 
to be the shape existing after inflation and after inplane edge loads 
are applied. All displacement quantities in the derivation are to be 
interpreted as being measured from this undeformed state. 

A rectangular xyz coordinate system is chosen so that the middle 
surface of the plate lies in the xy- plane as shown in figure l(a). In 
addition to the basic xyz coordinate system, two other coordinate systems 
are used in the analysis. These additional systems are designated x' 

y+, and x_, y_, z_ . The (+) subscript refers to the upper 

cover of the plate or the cover which lies on the positive- z side of the 
middle surface and the (-) subscript refers to the lower cover. This 
convention is used throughout the report. 

The primed coordinate systems are obtained from the xyz system by 
small rotations equal to the taper angles of the covers. Thus, the 
x + y + - plane is parallel to the upper cover and the x^y_[- plane is parallel 

to the lower cover. The coordinate transformations defining the primed 
systems are 


(i) 


as 

shown in figure 1. These quantities are assumed to be small ( i.e., 

b « 1, h « 1, and h h _ r « 1 ) and constant throughout the plate. 
* A > y j x ,y / 

DERIVATION OF DIFFERENTIAL EQUATIONS OF EQUILIBRIUM 
AND BOUNDARY CONDITIONS 


* = x + - z + — 


y = - z + 


x 

z = z| + x| -i- + y| 


h ,y ^ 


. x 


y = y_ + z, 


h V 
t y J 


h x h V 

z - z 1 - x ’ — y 1 — — 

" 2 2 J 


where the quantities h / 2 and h / 2 represent the taper angles 


Derivation of Stress- Displacement Relations 

It is desirable to derive the theory in terms of stress and dis- 
placement quantities which are associated with the overall plate con- 
figuration rather than local quantities associated with the individual 
covers because these overall quantities can be identified with familiar 
quantities in plate theory. The relationships between the local and 
overall quantities are discussed in this section. In addition, there 
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are derived relations between the overall plate stress resultants and 
the overall plate displacements. 


Strain-displacement relations .- Displacements in the individual 
covers of the inflatable plate in the prime! coordinate systems are denoted 
by u|, v+, and w+. These quantities are related to u+, v+, and w+, 

the displacements in the xyz coordinate system, by the transformations 



The strains in the individual covers are defined in the primed coordinate 
systems in the usual manner: 


€ x± 


Su± 1 /3w t 
Sx+ 2 \dx+ 






dwp 2 


Sv+ 1 

dy+ ^ W+ 


7 xy± 


Sll| dv+ ch'+ j 

dy+ Sx| S::+ dy+ J 


( 3 ) 


The following plate displacement quantities are now introduced: 


u 


u + + u 
2 


V 


V + + V 

2 


a 




w 


w + + w_ 


1 


e = 


v+ - w_ 


J 


(o 


The quantities u, v, and w are simply :he averages of the displace- 
ments in the upper and lower covers in the x-, y- , and z- directions, 
respectively. For small displacements, the quantities a and P are 


vo H vo tr 1 
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the angles between the drop cords and the z-axis in the xz- and yz-planes, 
respectively, or, in other words, a and p are the rotations of the 
drop cords during deformation. These displacements and rotations are. 
shown in figure 2 in their positive senses. The quantity e represents 
the unit deformation or "strain" through the thickness of the plate. 

Since the drop cords are assumed inextensional, e is not an independent 
quantity for small deformations but is related to a and p by the 
equation 



With the use of equations (l) to (k) the following sum- and- difference 
quantities necessary to the subsequent development may be calculated: 


x+ + € x- 

e x+ - e x- 

1 1 

y+ + € y- 

1 t 

e y+ " e y- 

' + j ' 

xy+ xy- 

7 ' - 7 1 

xy+ xy- 


In order to illustrate these calculations, consider the quantity 


€ x+ + e x~ 


8u|. l/^ w +\^ ^ U - 

8x_|. 2 \dx+ ) Sx_[ 2 y 8x J. / 


(6) 


obtained from the first two of equations (3). From equations (l) and ( 2 ) 


Su! 


8x^_ dx 


u 


.(x,y) + - 2 — w (x,y) 


dx + 8 


8xj 8y 


u + (x,y) + w + (x,y) 


Sy_ 

8x_! 


- lu + + 


w, 


( 7 ) 


Similarly, 


dw * 


1 x 


,y \ 


U V I 

2 + 2 + / 


( 8 ) 


and so forth. When equations (7) and (8) and their counterparts for the 
derivatives of u_[ and w^ are used and when it is remembered that 
h^ x and h y are constants, equation (6) becomes 
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h 


e x+ + e x- 


= u+,x + -p“ w +,x + u -,x 


12 


w -,x + -> w - 


h fx 2 


2 2 


r +,x + ^ v+,x. 


h >y ..2 u u .. „ , h ,x h ,y 

+ ^ v +,x " x w +,x u +,x _ t yU, x v +, x + g u +,x v +,x 


2 2 

2 h , x 2 h ,y 2 , 

+ w v + — — u v + — — v „ + h v w v u 

- j •*- 1). '> A 


4 -, x ,x -,x -,x 


+ h ,y w -,x v -,x + 2 


h ,x h ,y 


U -.x v -,x 


(9) 


The use of equations (4) then leads to 

h 


£ x+ 4 e x- = 2u ,x 4 ”fx 4 -r( he ,x + h ,x e ) 4 £( he ,x + h ,x e ) 2 4 “IT 


4 -f V 4 3J-( ha ,x 4 * x°f 4 #Hx 4 \J)‘ 


h ,x . x h ,x ; 


- ~2~ W ,x( h ,x a + ha ,x) - “H he ,x + h ,x e ) u ,x 


“ ~ 2 ^" w , x( h , x^ + h ^,x) " 2 ^( h€ ^ x + h ,x e ) v ,x 


+ - g - - u ,x v ,x + 'jj ,y ( h ,x a -• ha ,x)( h ,xP + h P,x) 


( 10 ) 


After substitution of equation ( 5 ) certain terms can be neglected in equa- 
tion (10). In the first place, terms which contain u, v, w, a, and p 
and/or their derivatives to higher than second degree are assumed to be 
small. Such terms ultimately lead to nonlinear equilibrium equations 
and are neglected in this theory. Second- degree terms in equation (10) 
lead to linear terms in the equilibrium equations and are retained here. 

Secondly, terms containing h^ x , h^ y , or h^ x h^ y are neglected because 
the taper is assumed to be small. 

When the approach just described is extended to the other sum-and- 
difference quantities the results are 


U) M VO f 
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e x+ + € x- 


= 2u^ x + wf x - |(h #x a x( a + w ,x) + h ,xP,xP + h ,y P ,x w ,x] 


€y + + £y_ = 2V^y + W^y - |^yP,y(P + W^y) + ^ y<» ytt + h f X tt , / W , £| 

7iy + + 7iy_ = 2u y + 2v x + 2v x w^y - |f, x a y (a + * x ) 



+ h ,yP,x(l 

3 + 

w ,y) + h ,xP,yP + 

h Y a 
y A y A j 


+ h v a v a 

} J 

+ h ,y p ,y w ,x] 


1 

- £ x- 

* ha , X + h ,x( 

a + 

w ,x) " h ,x w ,x u ,x 

- h v w 
y J y 

- 4 - 

- h|J ,y + h ,y( 

P + 

w ,y) " h >x w ,y u ,y 

- h v 
y j y 


- 7^. = ha y + h P>x + h^y (a + w >x ) + h >x (p + v y ) 

' h ,x( w ,x u ,y + w ,y u ,x) - h ,y( w ,x v ,y + w ,y v ,x) 


>( 11 ) 


Stress- strain relations - The strains in the covers are related to 
stresses in the covers through the orthotropic stress- strain relations 

N i± = A ll e x± + A 12 € y± 

^y± = A 21 6 x± + A 22 6 y± 

N xy± = A 55 7 xy± 

when the principal directions of the covers are alined with x+ and y ± 

axes. Note that the N' quantities, the stresses in the covers, are 
expressed in terms of force per unit length rather than the conventional 
force per unit area. If the covers are simply woven fabric, for instance, 
with the warp and fill alined with the and y' ± axes, respectively. 


■\ 


( 12 ) 
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then the coefficients of the strains in equations (l2) are related, to 
the usual orthotropic elastic constants for the cover materials by 


*11 = 


E W t 


1 - 


^WF^FW 


12 


PFW^ 1 


1 - 


^WF^FW 


"A 


*22 


^1 " 


Ejpt 

1 " ^WF^FW 
1 ' ^WF^FW 


Ajj — Gt 




(15) 


The Poisson's ratios are related by Pp^t = p^Ept; therefore, 
A 12 = A 21 • 


Stress- displacement relations .- It is appropriate, now, to intro- 
duce stress quantities N x± , N^, and N xy + associated with the 

xyz coordinate system and to define the following overall plate stress 
and moment resultants: 


X 

II 

N x + + N x _ 

M x = |(N X+ - 

Nx-) " 


II 

it 

V + N y- 

= I'V “ 

N y-) 

> (14) 

N xy : 

N xy+ + ^xy- 

**xy = rK N xy+ 

" N *y-L 



The assumption that the taper angles are small so that second-degree 
terns in h^ x and h^ y are negligible compared to unity justifies 


replacing 

equations 

N x± by N'±, Ny+ 
(l4) become 

by N^+, and 

K xy± 

by Ky±- 

Thus, 


N x + + N x . = N x 

N x + - 

V 

X- 

r 

II 



n ; + - n;_ = Ny 

N y+ * 

4- 

II 

> (15) 


N xy+ + N xy- = N xy 

V 

N xy- 

:>» 

i * 

II 

1 



t— I CT\ ( — I CT\ 
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The sums of the stress resultants in the upper and lower covers become 
plate stress resultants and the differences are associated with plate 
moment and twist resultants. These plate quantities are shown in fig- 
ure 2 in their positive senses. 

When the three pairs of equations (12) are added and subtracted 
there results 

N x+ + N x- = A ll( e x+ + e x-) + A 12 ( € y+ + e y-) 

N y+ + N y- = A 2l( € x+ + e x-) + A 22 ( e y+ + € y-) 

^xy+ + ^xy- = A 33( 7 xy+ + 7 xy-) 

N x+ " N x- = A ll( e x+ " e x-) + A 12 ( G y+ ” e y-) 

Ny + - Ny_ = A 21 (4 + - «x-) + A 22 ( e y+ " e y-- ) 

^xy+ " ^xy- = A 33( 7 xy+ “ 7 xy-) 

Equations (ll) and ( 15 ) are now substituted into equations (l6) and all 
nonlinear terms are neglected to obtain the following linear stress dis. 
placement equations: 
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Use of the Principle of Minimum Potential Energy 

The differential equations and boundary conditions are derived by 
using the principle of minimum potential energy. For the problem under 
consideration this principle may be stated as follows: When the plate 

is in equilibrium under external and body forces, the variation of the 
total potential energy with respect to variations in the five displace- 
ment quantities u, v, w, a, and 3 must be zero. The total poten- 
tial energy is the sum of the strain energy :n the covers and the poten- 
tial energy of the internal pressurized gas minus the work done by the 
external and body forces. The principle can be written 

6(nc + % - w) = o (18) 


where the three quantities 6IIc, SIIj, and 6W are expressed in terms 

of plate stress resultants and plate displacements. The variational 
operations indicated by the symbol 6 are performed in detail with 
respect to the displacements. When the coefficients of 5u, 6v, 5w, 

6a, and 63 are equated to zero there results a system of five linear 
partial differential equations and associated boundary conditions which 
govern the behavior of inflatable plates. 


Strain energy in covers .- The variation of the strain energy in the 
covers corresponding to arbitrary variations 6u, 6v, 5w, 6a, and 63 

of the displacements may be written as follows for a slightly tapered 
plate of arbitrary shape: 


611 c 


■JT\ 


K 


x+ oe x+ 


N ;* 6e y+ + 


N xy + ^W 


, d(xj,y+) 

d(x,y) 


+ 



t * 

+ 


y- 


+ N 


f * 

xy- 



dx dy 


(19) 


where 


1* i 0 1 

N x ± = N x -f + N^ 

t * ,0 . 

V " N y± + V 


__ » * » o , 

^xy± ” ^xy± + ^xy± 


>-l cr\ h o\ 
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or, since the Jacobians 
unity. 


^(x+jyt) 

d(x,y) 


of the transformations (l) are both 


&n c = JJ |(Nx+ + Nx!)s(ex+ + e x-) + |( N x+ “ N x-) 6 ( e x+ “ e x-) 


+ |(n y* + N y-) 6 ( e y+ + 4-) + |( N y* " R y-) 5 ^ e y + " 


+ i(Niy + + N xy-) 5 (yxy+ + ?xy-) + |( N xy+ " N xy-) 6 ( 7 xy+ " 7 xy-) 


dx dy 
( 20 ) 


The region of integration is the middle plane of the plate. Equa ion ( 
can be expressed in terms of plate stress and displacement quantities by 
using equations (ll) and the following equations analogous to equa- 

tions (15): 


» * » * , T * 
N x+ + N x _ = N x 


N 


» * 
'x+ 


N 


> * 


21VL 


- »• 


_ ^ 


- v* „ 


N 1 * + n’* = N* 

1N xy+ + w xy- xy 


N. 


» * 


2M- 


xy+ 


N. 


xy 


xy- 


When these substitutions are made, terms which consist of the product 
of inplane stress resultants (N*, N*, N* y ) and second-degree displace- 

ment quantities are retained because these stress resultants contain 
contributions which are independent of displacements (such as from 
internal pressure) and hence yield linear terms in the differential 
equations. Elsewhere, second-degree displacement quantities are neg- 
lected. With the use of equations (ll) the expression for the variation 
of the strain energy in the covers can be written as follows wi s lg 
rearrangement : 
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6JI f 


= 11 


(n*8 

r~ 

1 2 

u „ + — v „ - 

hh x 

— 7— ■ — a (a + w \ 

hh,x 

y a p 

hh v "1 

l x 

,x 2 

4 ,xV a + w ,sJ 

4 

"T“ 


+ N* y 5 


hh 


u 


- - 

hh ,y 

■ V 21 3 

b > 

(b + 

, hh 

w ) - ,y 

a ,y a “ 

hh v 

) X 

\ ,y 

yv p 

>y) 4 

4 ,y 

v + 

i 

>> 

*N 

X 

> 

hh v 

a ,y( a + w ,x] 

c 

c 

>> 


4 

1 ‘ 4 


hh Y 

y A 

aw 
>x , y 

^y 

X 

P 

i 

D 

+ M 5a 
x y : 

4 

4 

p W 

>y ,x 


y 


r— 0 0 - 

+ M y 6f3 ,y + ^y 5 ( a , y + p ,x) + s x 6 ( a + w , x ) + s y 5 (0 + ^yjj^ dy ( 21 ) 

where S x and. Sy are defined by 

s x = ^-( h , X M X + h ,y M xy) 

Sy = h^x^y + h ,y M y) 


( 22 ) 


These quantities represent that portion of the transverse shear which is 
carried in the covers by virtue of the taper. 


Potential energy of internal pressure .- Phe variation of the poten- 
tial energy of the pressurized gas inside the inflatable plate is the 
product of the magnitude of the internal pressure and the change in 
volume due to virtual displacements 5u, 6v, 8w, 6a, and 60. The 

increase in volume of the plate due to any set of small displacements 
u ±> v ±> an< i w+ in the covers and along the edge, denoted AV, is 

calculated from 


AV = - 



CU,Tl)d| dT] - 


T 


h dx dy 


(23) 


In the first integral r\ } and £ are rectangular Cartesian coordi- 

nates locating the final positions of points on the surface of the plate 
and are functions of certain of the displacements u+, v+, and w+. 

The surface integration is performed over the entire outer surface of 


hi C\H C\ 
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the plate (over the upper and lower covers and along the edge) in such 
a manner that the surface is oriented with positive inward normal. (See 
ref. 1.) The second integral is performed over the middle surface of 
the plate at its initial position and represents the volume before defor- 
mation. The details of this calculation are presented in appendix A. 


The variation of the potential energy of the internal gas can be 
written BJI-j- = -p5 AV and the results of the calculation in appendix A 

yield 


&H-J- = -pB 


u 


+ v 


+ u 


* ,y 


- u 


+ x p 

12 N > X > 





(2k) 


In this expression quantities which are higher than second degree in u, 
v, w, a, and (3 or their derivatives are neglected. The double inte- 
gration is performed over the unde formed middle surface of the plate. 


Work of external and body forces .- The plate may be loaded by dis- 
tributed forces, inertia forces, and forces and moments acting on the 
boundaries. The variation of the work of these forces may be written 
as 


5W 


-M 


q^u + lyBv + q. 


u2 .. \ 

,5w - pu5u - pv6v - pw5w - pa5cx - -jj- ppbpjdx dy 


+ 




N Nx 5u 


Njjy6 


v + V n 5w + 


M^yba — ds 


(25) 


As before, the double integral is performed over the middle surface of 
the undeformed plate. The line integral is around the boundary of the 
undeformed middle surface and positive in the clockwise sense looking 
in the positive z-direction. (See fig. l(a).) 


Equilibrium Equations and Boundary Conditions 


The previously derived parts of the variation of the potential 
energy are now collected, and where required variational operations are 
performed in detail. Equation (l8) thus becomes. 
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sn = 5n c + &Hj - &w = 



n: 


X 


Bu 


hh 


mi 

>x + w , X 5 w X - -jp- a x (6a + &v x ) 


hh ,x, v hh hh x 

— (“ + W ,x) 5 V - ~t~ P ,x« - -f 2 - TO, X - 


hh 

~T — 3 x Bw . 
4 > A > - 


hh y 

> J ctn 

— w ,x 5 P,x 


+ N- 


y 


5v 


hh 


,y + w ,y 6w ,y ' -j~- 3,y(&P + 6v >y ) 


^(3 . » ,)*„ - iz a y Sa - ^ ^ ^ „ y6 „ 


4 

hh 


► y ,y 


, x 

— w 5a 

4 *y >y 


+ N 


*y 


5ii + Bv + w bw + w , r 5v 


,y 

hh 


,x * ",x— ,y ' w ,r,x 

hh 


hh^ x hh x hh 

— a y (8a + 6w^ x ) - -j^(a + v JSa y - -^1 0^( B p + 6v^ y ) 


4 

hh 


4 

hh 


^* v)6(jI .^ v .!V W)r !V v6l 


, A 

w &a 


4 >y > x 4 

hh 

IT v .^,y 


hh y hh 

-a> ll a6a 


y " T~ > x ,y 

hh 

TT p »y 8w .x 


“x^.x + V#,y + *«xyH y + S0, x ) 


6u + 6v + u y 5v , 

>* >y > x 


+ S x (5a + 5w^ x ) + S y v &3 + 5w^ y ) - ph 
+ v ,y 6u ,x - u , y 8v ,x - V 5u ,y + lI(V 6f3 ,y + 0,/v - 

P,x^ a , y ) - a ^ a ~ 3®3 - aBw^x - w^ x Sa - 05w y - w y &pj 


q x B u - Q. y Sv - q z 5w + pu5u + pv6v + pw8v + ^ pa5a + ^ ppsp^dx d y 
4 (%x&u + % y &v + V N 6w + MjjySa - M Nx 83^ds = 0 


(26) 


>A 0 \ rH Q\ 



Q 


17 


The terms with partial derivatives of varied quantities are integrated 
by using Gauss' theorem which can be written as follows (see ref. 1): 


K* + ^ ' - *fi ( f i* + * 


dj ' ’*• J \ 6n 13 8n J 

The operator 6/8 n denotes differentiation in the direction of the 
inward normal. For example, a term like Jj ° N*5u x dx dy in equa- 
tion (26) is handled as follows: 


( 27 ) 


r 


( N x 5u ),x dx dy = - y N*8u ds 


or 


N x 8u ,x dx dy = - O N*5u ^ ds - JJ N*^ x 5u dx dy 


(28) 


( 29 ) 


When Gauss' theorem is applied in this fashion to terms in the double 
integral of equation ( 26 ), there results 

“ JJ ^( N x,x + N xy,y - Ph, x P h ,x v ,y + P h ,y v ,x + *x - P^Bu + ( N y,y 

+ N* y ^ x - Ph >y - Ph, y u x + Ph >x u^ y + q y - pv)6v + |(n* >x + N* y>y )w ;X 

+ ( N yj,y + N xy,x) w ,y + N x w ,xx + ^xy^xy + N y w ,yy + S x,x + s y,y + P h ,x a 

hh x C. 

+ Ph^ y P + Ph(a^ x + P^ y ) ^ ( N x,x + N xy,y) a ,x + ( N y,y + N xy,x) a ,y 


+ N x a ,xx + 2N xy a , xy + N y a ,yy 


hh 


( N y,y + W xy,x)P,y + ( N x,: 


+ ^xy,y)P,x + N xP,xx + ^xyP,xy + ^yP, 


yy 


+ q z - pUj-8w + <jM^ x + M*. 


<y,y 


hh 


s x - ph(a + W jX ) n“( h ,xP,y - h jy 3^ x ) 7 ^— ( N x,x + N xy,y)( a + w ,x) 


(Equation continued on next page) 
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*| till 

+ ( N y,y + N xy, x) w , y + N x w , xx + ^xy^xy + N ;' w ,yyj i~( N i 


h 2 •• 


_ __ p a^5a + <jMy jy + M xy ^ x - S y - ph (f3 + y ) — (^y^x “ h ,x a ,y) 

i^~|( N y,y + N xy,x)(p + w ,y) + ( N x,x + N xy,y) w ,x + N x w ,xx + ^xyW, j 
+ N y w ,yyj 4 ( N x,x + N xy,y)P ■ ~ P^Bpjax dy 


Nx rn + N xy In " ph § " phv .y S? + ph V % + %x) 5u 


, (n* dy * dx , dy , dy dx - \ 

\ y dN xy dN ‘ P dN ■ P > x dN + P > y dN + N Nyj 5v 

+ few x £ + N^v y |y + N xy w x & + N* ■ V y ^ + S x |E + 8 y & 

[_ x ,x dN y > y dN xy ' x dN xy > y dN x dN ^ dN 

+ Pta | + PM | ♦ h - ^(4«, x 1 ♦ | * »V,y £ 

+ "W 1) - | * »V,x | * »V,y 1)1- 


f r^x Sy — hh 

+ r* dN + dN + %y ‘ T 


) x M * / \ d> * dy * dx 

— N x (a + v v )— - + N^w — + N , r w — 

L xV > x 'd]v y ,y dN xy > y dN 


* ”,x)|] - 


+ B * a 2i\ . £S!f p 

dN x y dN/ 12 \ P > y dN 


Sy \U„ . L dy . „ dx w ... N ^y . dx 


- »,« i - - 7 ^ ♦ “,y)i ♦ &,* i 


n *y“,x £ + »xy( p + w ,y)f|] - ("> £ * "V 

ir( a ,x | - a .y |)W ■ ° 


CTn 1 — I ON 
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In writing equation ( 50 ) second- degree terms in h x and h y are 

neglected, and use has been made of the fact that h Y and h __ are con- 

y A y y 

stants. By the usual arguments of the calculus of variations, the coeffi- 
cients of the quantities 5u, 6v, 6w, 5a, and 5(3 in the double integral 

can be equated to zero individually. This procedure gives a system of 
five partial differential equations of equilibrium- The first two equa- 
tions describe the equilibrium conditions for the inplane forces, the 
third is the lateral- force equilibrium equation, and the final two are 
moment equilibrium equations. 


The last three partial differential equations, those associated with 
lateral forces and moments, can be simplified by some additional manipu- 
lation. The lateral- force equation is simplified by using the two inplane 
force equilibrium equations to substitute for the quantities N* + N* __ 

-K- * x y y 

and Ny^y + N X y^ x . Terms in the resulting expression are neglected if 

they contain second-degree quantities in h Y and h „ or nonlinear 

y ^ y y 

quantities in the displacements or their derivatives. In addition, terms 
which contain products of an external loading and a displacement quantity 
are nonlinear and are dropped. The moment equilibrium equations are 
simplified by using the inplane force equations and also by substituting 
for the quantity N^w^^ + 2N X y w , xy + ^y v , yy ? rom the simplified lateral- 

force equilibrium equation. The same sorts of quantities are neglected 
as mentioned previously. In order to neglect all such quantities, the 
assumption must be made that the inplane stress resultants (n*, N*, N*y) 

are not of a higher order of magnitude than the quantity ph. The 
simplified equations of equilibrium which result from these manipula- 
tions are 


* 4 " N — 

*y,y 

P h ,x 4 + \y) 

+ P h ,y v ,x + 

= pu 

(31a) 

+ w* 

+ n xy,x 

ph , y 4 + u ,x) 

+ ph x u y + q y 

= pv 

(311) 


4 " £( h ,X a + h ,y P )] >xx + N y(^ ' 4 X<X + 

+ 2Nxy|w - x a + h^y3)J + S x>x + Sy >y + ph >x (a + w jX ) 


+ ph^ y (p + v^ y ) + ph(a x + 3 >y ) + q z = pw 


( 31 c) 
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M y y + M v 
v ■ 


- S x - ph 


a + w ,x - r ( h , x a + h ,y3), : 


hh Y ,2 

+ -j^-(q z - P«) = \ 


(31d) 


*V,y + ^y^x s y ” ph 


+ M , y - ^( h ,x a+ h ,/0h 


hh __ p 

+ ~ r~(q z - pw) = V pP 


4 

( 31 e) 


Note that in equation (31c) the stress resultants N*, N*, and N* y 

be replaced by N°, Ny, and N^y in order :o linearize the equation. 

r 

h 


must 


Nonlinear terms like N 


w - h( h ,x a + h ,yP) 


are assumed to be small 


, xx 


compared with the linear terms and are neglected. 

The appropriate boundary conditions are obtained by equating to 
zero the line integral in equation (30). The usual arguments of the 
calculus of variations lead to the following boundary conditions on the 
boundary of the plate middle surface: 

u = 0 "I 


or 


Nx - ph(l + V } y) cos(N,x) - (Nxy + ph x ) COS (N, y ) = % x 


( 32 a) 


or 


v = 0 


n: 


P h (l + u^ x ) 


cos(N,y) - (N* y + pfc cos(N,x) = % y 


( 32 b) 


w = 0 
or 


-fekx ' jj( h ,x a + h ,yP) 
+ phalcos(N, x) - 


+ N 


xy 


w ,y - r(fx a+ h ,y3) 


,y 




+ 


J 


V - £' h 


',x a ^ h ,/) y ] 


) ( 32 c) 


+ N 


xy 


W ,X - £( h ,x a + h ,yP) 


+ S, 


+ php|co£ (N,y) = V N 


J 


vo h vo Ir* 
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P^l R - £ 
12 p >y u 


a = 0 
or 

- 

ph^ . h 

+ 12 ~ P ,x " 4 


(3=0 

or 


M . a h 
W 12 > x 4 


■N 


^x^,x( a + w , x) + ^xyO 1 , y a + x w , y)|j cos (N> x ) - 

, 'xy h ,x( a + w ,x) * K ?( h ,y a + h ,x",ym cos (“-y) = «Ny 


J 

( 32 d) 


N^x(P + w ,y) + N xy( h ,xP + h , y w , x) |c°s( N ,y) + ^ l 


V 


ph^ h 

+ ~ “ y " k 


-is 


N xy^, y(P + w ; y) + N x( h ,x p + h ,y w ,x) f cos ( N,x ) M ] 


V 


Nx 


J 

(32e) 


In these equations the relations 


^ = cos(N.x) 

SN ' 


^ = cos(N,y) 




(33) 


are used where (N,x) represents the angle between the inward normal and 
the positive x-axis and (N , y) represents the angle between the inward 
normal and the positive y-axis. The differential equations (31) and 
boundary conditions (32) govern the behavior of inflatable plates. In 
addition to these equations the linearized stress displacement rela- 
tions, equations ( 17 ), are required to complete the system. 


Limitation of Linear Theory 

The range of validity of linear theory certainly depends on the 
specific problem under consideration. In some problems linear theory 
might be valid up to substantially larger deflections than in other 
problems. In order to obtain a rough indication of what this range of 
validity is, one simple problem is investigated in appendix B. The 
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lateral deflection of a long plate simply supported on its long edges and 
under uniform load is considered where the simply supported edges of the 
plate are prevented from moving together. Calculations made by a simpli- 
fied nonlinear theory are compared with calculations made by the linear 
theory. The results show that the linear result is within 5 percent of 
the nonlinear result if the following condition holds: 

2 

w 

(J + l)-i < 0.05 


where 



where wp is the lateral deflection at the center of the plate, and 

where L is the distance between the simply supported edges. Of course, 
this result is strictly applicable to this one problem only but it does 
give a rough indication of what might happen in other problems. 


SOME SOLUTIONS FOR RECTANGULAR PLATES 


Solutions are now presented to some specific problems in the deflec- 
tion and vibration of constant- depth rectangular inflatable plates sup- 
ported on all edges. Consideration is given to static deflection under 
uniformly distributed lateral load and vibration modes and frequencies 
of both simply supported and clamped plates. A comparison between 
results calculated from these solutions and some experiments conducted 
on a square inflatable plate is presented ir reference 2. 


Consider a rectangular plate of constant depth supported along 
edges x = 0,a and y = 0,b. For the case where no edge loads are 
prescribed, the boundary conditions, equations ( 52 ), are as follows: 


on x = 0, a. 


or 


u = 0 

< - pH 1 + \y) * 0 



v — 0 


N*y + phU y = 0 



(5^a) 


(5^b) 


or 


1-3 ON H < 7 \ 
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or 


or 


or 


and on y = 0,b, 


or 


or 


or 


or 


or 


w = 0 


N x w x + -^xy w , y + Ph a - 0 


a = 0 


M * - h- 


= 0 


) 




0 = 0 


«y- 


ph^ 

12 ~ 


a 



(54c) 


(54d) 


(54e) 


(55a) 


(55b) 


(55c) 


(55d) 


(55e) 



2b 


If the u and v accelerations are neglected and only lateral loads 
are considered, the differential equations I 3l) become 


N x,x + N xy,y “ 0 
N y,y + N xy,x = 0 


N x w ,xx + N y w ,yy + ^xyVy + Ph(a x + P, y ) + <l z = pv l {%) 

Mx.x + Mxy.y " P h ( a + u ,x) = X ° S 
»y,y + M xy,x - PKl 3 + “,y) " T pg 

The portions of the stress resultants which are independent of dis- 
placements must satisfy the following differential equations: 

N° + N° =0 
iM x,x + "xy,y u 

N° + N° =0 

y,y xy,x 

and, if u and v are not prescribed on the boundaries, the following 
boundary conditions must be satisfied: 


on x = 0, a. 


and on y = 0,b, 


- ph = 0 


N° =0 
xy 


N° = 0 


- ph = 0 


The solution to this system is 


N° = N° = ph ^ 


N° =0 
xy 


*— 1 CTn i — I 
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The first two of equations (56) contain only u and v and the remaining 
three contain only w, a, and p. Since for the particular problem under 
consideration the lateral deflection is desired, the fact that the u 
and v equations uncouple from the others means that it is necessary to 
consider only the last three of the differential equations (56) and the 
appropriate boundary conditions. These differential equations reduce to 


ph(a + w^)^ + ph(0 + + q z = pw ^ 

Wx,x + ^xy,y ' P h ( a + w ,x) = P ' d / 


2 .. 


M- 


+ M xy,x - ph(p + v y ) = — pp 


and the boundary conditions become: 
on x = 0,a, 

w = 0 


or 


ph(a + WjX ) = 0 


or 



(58) 


(59a) 


( 59 b) 


or 


and on y = 0, b, 


or 


P * 0 


M- 


xy 




w = 0 


ph (3 + w ,y) = 0 



(59c) 


( 40 a) 
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or 


or 


a = 0 


M- 


•xy 




3=0 


ph 


Mv - — a Y = 0 

y 12 ,x 


(40b) 


(40c) 


Finally, from the last three of equations (17) the moments are written 
In terms of the rotations 


= ^( A ll a ,x + A 12 P ,;v) 

*V = + A 22 P ,;/) ) 

= ~ A 33 C a , y + P ,x) j 


(41) 


The differential equations ( 58 ) have the same form as the well- 
known equations for the lateral motion of a uniform plate with trans- 
verse shear flexibility and rotary inertia. (See, for example, ref. 3 .) 

In this case the quantity ph plays the role of the transverse shear 
stiffness per unit width. The salient conclusion of this investigation, 
therefore, is that an inflatable plate can be considered as a particular 
type of sandwich plate where the pressurized gas acts as the core material. 
Note, however, that certain terms peculiar to inflatable plates appear in 
the boundary conditions (39) and (^O). These terms are associated with 
the deformations of the edge walls in the presence of internal pressure. 
Consideration should be given to the important of these terms before 
sandwich- plate solutions are applied to inflatable-plate problems. 

When equations (Ul) are substituted into the differential equa- 
tions ( 38 ) there results 


v-t CTn 1— I G\ 
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P h ( a + w ,x)^ x + P h (P + w ,yh 


+ 4- 


pw "I 


h 2 *,2 0 h 2 

2 ( A ll a ,xx + A 12P,xy) + 2 Mtyy + P ,xy) " P h ( a + w ,x) = ~ 
^ 2 * ( A 21 a , xy + A 223,yy) + %• A 33( a ,yx + P ,xx) ' P h (P + w ,y) = £ ^“ 


a 

P 


(42) 


The differential equations are now expressed entirely in terms of the 
plate lateral displacements w and plate rotations a and p. 


Simply Supported Rectangular Plates 


For the case where all boundaries of the plate x = 0, a and y = 0, b 
are simply supported, a reasonable set of boundary conditions can be 
written from equations (39) and (40), with the help of equations (4l), 
in the form 


and 


w(0,y) = w(a,y) = 0 
p(0,y) = P(a,y) =0 > 

a x (°>y) = a = 0 

j A J 

w(x,0) - w(x,b) = 0 
a(x, 0) = a(x,b) = 0 / 

P (x, 0) = p (x,b) = 0 

y J > J J 


(43) 


(44) 


Static deflection under uniform lateral load .- The response of the 
plate to a uniform lateral load is governed by equations (4-2) with 
terms w, a, and p set equal to zero and the term q z regarded as 

constant. It is seen that the boundary conditions (43) &nd (44) are 
satisfied by displacements and rotations of the form 
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m=l n=l 
00 00 
•■IE 

m=l n=l 

co oo 

'•T.Z 

m =1 n =1 



(^5) 


When equations (4-5) are substituted, into equations (4-2) reduced as just 
indicated, there results 





Equating to zero the coefficients of the Foirier series in the last two 
of equations (46) yields equations for and in terms of v mry . 

The constant term q z /ph in the first of equations (46) can be expanded 

in a double sine series in the interval 0 = x ^ a, 0 ^ y ^ b to yield 
the equation 
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w 


2i + 

£) 

L Ann 

m 

U 

A 

T “ 

Jt 

a 

(— + 


Ct mn 

-i- " ■ 

m 

u 

b 2 J 

n 

a 


fimn n _ l6 _1_ 
jt b ^ ph mn 


Ulii -LJ. _ 

Tb " 



(^ 7 ) 


When the expressions for and. ( 3 ^ in terms of are substi- 

tuted into equation ( 47 ), the resulting equation can be solved for v^. 
This procedure leads to the following expression for in dimension- 

less form (for m and n odd): 


w. 


mn 

h 


Rq 


mn 1 - 


^") a mn - 77^- (fnm + (^ 2 ~ n 9% 5 _ 

Kin / SAnn ^ 


m 


where 


* = b 
^ a 


a = 




* 2 p 


R = 


4 b 2 

n 2 h 2 


P = 


Ph 

A 11 


Ann 4 


( R P + b nm)( R P + c ran) “ AL 


,2 2 2 

Ann = A m + n 

Ann = 2(x 2 m 2 + A 55 n 2 ) 
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and. where 


c mn = 2 ( A 22 n2 + A ^y^m 2 ) 

V = 2 ^ nn ( A i2 + A 5j) 

222 / — — \ 
e mn = A m n (1 + A 22 - 2A 12 J 


a 12 


A 12 

A 11 


A 22 


A 22 

A 11 



= ^33 
' All 


For either m or n even, all coefficients 


w mn> ot rnn-' ^ 


Pmn vanish. 


Natural frequencies of vibration .- The free vibration of a s imp ly 
supported rectangular plate is governed by differential equations (42) 
with q z set equal to zero and boundary conditions (eqs. (43) an d (44)) 

assumed to hold for all time. The functions (for m = 1, 2, . . . and 
n = 1, 2, . . .) 


w = 


a = 


P = 




iayt . mux . nny 
w mn e sin -r- sin ~rr 


ia>t mnx . mty \ 
W cos — Sin ) 


lent mnx mity 


3mn e sin — ccs 


b J 


( 49 ) 


are seen to satisfy the boundary conditions, equations (43) and (44). 
When equations (49) are substituted into equations (42) (with q z = 0) 

there results (for m = 1, 2, . . . and n = 1, 2, . . . ) 


G\ H Q\ 
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or in dimensionless form 


* 

l[p(AV + n 2 ) - k 2 ) 

pAm 

pn 


^mn 

h 

- 

pAm 

pb hn 2 f 2 2 - 2 k 2> \ 

~h + 2Tl A m + A 33 n * TJ 

|2- >mn(A 12 + Ajj) 


a tnn 

it 


pn 

lb - ^(*12 + A 55 ) 

't * ♦ *3!^ • 4 ) 


Pmn 

n 


? p 

2 pen 

where k = is the frequency parameter. When the determinant of 

An* 2 

the coefficients of w^, P mn I s equated to zero and after 

some algebraic manipulation, the following equation results: 


k^ + 


_ (2R + Qmn) + ^mn + Cmn 

+ ( b mn c mn ' d mn) k ~ P 2Rf mn " P®mn ( b mn c mn " d mn) = 0 
where the only symbol not previously defined is 


|p R ( R + a mn) + p( r + a mn)(l ) mn + c mn) 

( 52 ) 


f = A^b.. 


.jnn - » - -mn - » + 2Amnd mn 


+ n c„ 
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The solutions to equations (52) define the natural vibration frequencies 
of the simply supported plate corresponding to the mode shapes given by 
equations (49). 


Clamped Rectangular Plates 

Suppose, now, that instead of boundary conditions as given by equa- 
tions (4j) and (44), the following boundary conditions are prescribed 
for a rectangular plate: 


and 


w(0,y) = w(a, y) 
a(0, y) = a( a, y ) 
P(0,y) = P ( a, y ) 



w(x,0) = w(x,b) = 0 

a(x,0) = a(x,b) = 0 

P(x,0) = P(x,b) = 0 


(53) 


(54) 


These boundary conditions are appropriate for a plate with all edges 
clamped. The differential equations (42) must now be solved in conjunc- 
tion with conditions (53) and (54). 


Static deflection under uniform lateral load .- In equations (42) 
set w, a, and p equal to zero and regard q z as constant. The fol- 
lowing deflection and rotation functions satisfy the boundary condi- 
tions (53) (54): 


w 


-Z Z “■ 


A 


. mjtx n Tty 


rnn sin s:.n 

a b 


m=l n=l 


a -Z Z 


rrntx nny 

Omn sin — s: - n IT 


m=l n=l 


E r* . mnx nity 

2^ Pmn " — 


m=l n=l 


( 55 ) 


vo h vo tr 1 
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that is, w, a, and. (3 are zero on all boundaries. Substituting equa- 
tions (55) in the reduced differential equations (42) yields 


00 00 


r 

. ~ w mn 

(“f * ("f 

l 

— 


mrtx 

a 


nrty 

b 


sin sin + cy ™ cos = sin 


mrtx 

a 


nrty 

b 


, n ELJT . mrtx 

+ Pmn — sin — b 


nrtjrV 


ph 


E V n } , mrt mrtx , nity 

L < - Wh — cos — sin -JL - ^ 


m=l n=l 


ph + A 


2 2 


l H — + A 33 — )- 


n 2 * 2 \h 2 


. mrtx . nrty 

sin sin — - 

a b 


> (56) 


0 

, n mrt nrt h { /. A \ mrtx nityl _ 

+ Pmn — — — ( a 12 + AjjJcos — cos = 0 


00 00 


mitx nrty 

r*r\ c - V 

b 


m=l n=l 


V V* ) „ njt mjtx ^y mrt nrt h 2 / A . \ mrtx 

1 L \ ^ ~ ~ T + T T T( A i2 + a 33) cos — cos 

/. n 2 « 2 

r 22 " 


~ 3 m n 


ph + 


m 2 rt 2 \h 2 


33 T 2 T 2 


. mrtx . nrty V 

sin sin — - > = 0 

a b f 


The Galerkin method is now applied by multiplying through all of equa- 
tions ( 56 ) by sin sin 4 ^ (where i and j are integers) and 

a b 

integrating over x from 0 to a and over y from 0 to b. The 
resulting equations are written in dimensionless form as follows: 


v i j fh\ 2 (.2.. 2 2\ n [(-D* ~ - l] 

—lb)" ^ X + J ) + 1 mA Vim + Z np in r Jn = ~ 2 TJ 

m=l n=l ph,< 


Z mX IT 7im ' (§ + S“( l2x2 + A 33-3 2 ) °ij + ^( A 12 + a 33 )Z Z mXn, W'im7jn = 0 

m=l n=l 


)( 57 ) 


. y , CO 00 rn p . “1 

5 Z " T >> * i( S « * *53)2 £ ■‘V* - f * §-(%/ * S»A 2 ) ■ 0 

n=X m=l n=l L J 



5 ^ 


where (using 


r and. s 


as integers) 


(-l) r - s - 1 (-l) r+S - 1 

y = + 

'rs r - s i + s 


(r f s) 


7 rs =0 (r = s) 

Examination of equations ( 57 ) shows that the unknowns w^j are nonzero 
only for odd values of i and j, the unkncwns <Xj_j are nonzero only 
for even values of i and odd values of j, and the unknowns ( 3 ^ j are 

nonzero only for odd values of i and even values of j . For a first 
approximation, truncate the system (57) ty taking i = 1, 2 and j = 1, 2 
and summing over m = 1 , 2 and n = 1 , 2 . Equations ( 57 ) reduce to 



and three additional sets of equations which are homogeneous by virtue 
of the symmetry of the loading and, hence, yield the trivial results 
w^^ = W21 = w^^ = OLr^i “ 0*12 = CL22 = Pqq = Pop = = ^ * The last two 

of equations (58) can be solved for 021 and ( 3^2 in terms of wpp/h. 

Substituting these expressions into the first of equations (58) and 
solving for w-q/* 1 results in 



Rq 


1 - 


16 

9 n 


R 2 ? 2 
2 A 


a ll 


16 Rp L 

9 * 2 A V 


,2 

A c 2 + 


, 52 

b 21 

9« 2 


Ad. 


22 , 


( 59 ) 


where 


A = 


1 

4 


(Rp + b 2 l)(^P 


+ C 12 ) •• 



More accurate approximations can be obtained by taking more terms in the 
series for w, a, and [3, but the number of simultaneous equations which 
must be solved increases rapidly. 
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Natural frequencies of vibration .- For the free vibrations of a 
clamped rectangular plate assume deflection and rotation functions as 
follows : 


00 00 


w - L T, 


icot . mnx . mty 
mn e sin — sin 


m=l n=l 


CO 00 

V V - icot . m«x . mty V 

a * L L v 5in — sin b > 

m=l n=l 


p = 


m=l n-1 


icot . mnx . mty 


Pmn e sin 


sin 


J 


(60) 


When equations (60) are substituted into differential equations (42) with 
q z set equal to zero there results 


1 1 \ ph ' \ a 2 + t 

m=l n=l 


, mitx . n*y 
w mn sin — sin — 


+ ^ = cos SS sin 2S + - 3mn M sin 5S cos 22\ = 0 


OO 00 


m=l n=l V_ 

■> 

t mrt n« h 2 / A » \ mitx mryV 0 

+ Pmn - V -( A 12 + A 33) COS IT C0S T “ 0 


ph 2 2 h h 2 (, m 2 it 2 . n 2 fl 2 \ 
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Application of the Galerkin method as used in the previous section leads 
to the following system of equations in dimensionless form: 


Ar 2 

b 



00 00 "v 

p 2_, “ P Y n " P - 7 on = 0 

m=l n=l 
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00 00 

+ ^(^12 + Y mM ^ 7im7 j n = 0 

m=l n=l 


} (62) 


- V* ln h It 

P L n IT 7 J n 2h \12 + 

n=l 


55, 


)E Z “K.v im 7 jn 


m=l n=l 


+ 
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Limit consideration to the first four terms of each of the series 
in equations (60); thus, the integers i anl j in equations (62) are 
taken to be 1 and 2 and the summations are carried out over m = 1, 2 
and n = 1, 2. It is found that only certain terms couple. The w-q 

coefficient couples with and p-^, ^qn couples with a gg and 

$H, w^ couples with 34^ and and finally w^ 2 couples with 

<xi2 and $21* The frequencies for these various modes within the 
framework of the 1, 2 approximation are given by the following equations: 


For the w-^ mode, 


2^4 (k 2 - p(a 2 + ljj 
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' 3 P 


- ^ pA 


8 - 

3 P 


k 2 - — - ^-^(4A 2 + A,,) 
4b h 2b ' 33/ 
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Q b V 12 33/ 


9 1 


^ l *( A 12 + a 55 ) TT\~ k2 ' T " IT ^55 + UA 22 ) 


= o (65) 
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for the w-^2 mode, 
^[k 2 - i(A 2 + 4 


pA 


8 - 
5 P 


-|PA 


! p 


k 2 - - 25 l( 4 A 2 + UA,,) 

4b h 2b ' 33/ 

-f | X ( A 12 + ^ 35 ) 


f | A ( A 12 + A 55 ) 

& * 2 - ! - #^5 ♦ **) 


= 0 (64) 


for the V2]_ mode, 

- p(4X 2 + l3 
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- T - It( a2 + ^ 3 ) ?S*&2 + S ») 
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and for the mode, 


!^[ k 2 . p(4A 2 + 4)| 

f'pA 

8 - 

3 P 

f 5A 

^ k 2 - ^ - *^(x 2 + 4 a,,) 
5b h 2 b ' 33 / 

■ 1 4 s i2 i *55) 

i p 

* 1 a ( a 12 + a 33) 

& * 2 - $ - 


= 0 (66) 


In each case the smallest real positive value of k which makes the 
determinant vanish corresponds to the desired natural frequency. 


TOTAL POTENTIAL ENERGY EXPRESSION 


The derivation of the theory has been carried out using the poten- 
tial energy in varied form. It is of interest to write the total poten- 
tial energy expression because such an expression is useful in obtaining 
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approximate solutions. The potential energy is assumed, to he measured 
from the condition of the plate after it has been inflated and after 
inplane edge loads % x and % y have been applied. In this state 

the displacements and rotations u, v, w, a, and (3 are assumed to 
be zero and the quantities N°, N y , and N^ y are assumed to be known 

and independent of u, v, w, <x, and |3. Then, additional loads q x , 
q y , q z , V N , M Nx , and % y are applied. For this situation, the total 

potential energy for the small deflections of a tapered inflatable plate of 
arbitrary planform shape can be written in terms of displacements and 
rotations as 



- P h ^,x + v ,y + u ,x v ,y - W ----- aw , 

+ l|(“,xU,y - “,ytJj - ' l A** V 

- ^ (% x u + % y V + V N W + %ya. - % x p)ds 


( 67 ) 
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L 
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In writing equation (67) terms of higher than second degree in displace- 
ments and rotations or their derivatives have been neglected. In addi- 
tion, inertia effects are neglected. Variation of equation (67) with 
respect to u, v, w, a, and p yields with a little manipulation the 
differential equations (31) (with the inertia terms set equal to zero) 
and the boundary conditions (32). 

In many problems only lateral deflections are desired and inplane 
displacements are not of primary concern. Then equation (67) can be 
simplified for use in obtaining approximate solutions by dropping all 
terms containing u or v or their derivatives. 


For a rectangular plate with edges x = 0 , a and y = 0 , b the 
double integration in equation (67) is, of course, performed over x 
from 0 to a_ and over ^y from 0 to b. Ety noting that on the 

edge x = a, N Nx - N x , % y ^= N xy , V N = V x , M Ny = M x , and M Nx = -M, 
on the_edge x = 0 , H Nx = -N x , N Ny = -N xy , V N = -V x , M Ny = -M x , 

%x = ^xy> 8X1(1 so forth, the line integral in equation (67) becomes 





1+0 


the double integration in equation (67) can be performed first over y 

from 0 to 
gral becomes 


from 0 to c(l - and then over x fron 0 to Z. The line inte- 


%x u + %y v + V N W + M Ny a " M Nx^) ds 
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dx 
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n l ^ _ _ 
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CONCLUDING REMARKS! 


A linear theory has been developed which describes the behavior of 
inflatable plates such as Airmat. A rough indication of the range of 
applicability of the linear theory is obtained in an appendix by carrying 
out a nonlinear analysis of a simple problem and comparing the results 
with the linear solution. The theory turns out to be essentially the 
same as sandwich-plate theory in which transverse shear deformations 
are taken into account. The internal pressure in the inflatable plate 
is analogous to the transverse shear modulus in the sandwich plate. 

Some simple static lateral deflection and vibration problems are solved 
for rectangular inflatable plates of constant depth. 


Langley Research Center, 

National Aeronautics and Space Administration, 
Langley Field, Va., June 21, 1961. 
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APPENDIX A 

CHANGE IN VOLUME OF INFLATABLE PLATE 


Analytical Approach 

In order to calculate the work done against the internal pressure 
In an inflatable plate it is necessary to know how the volume of the plate 
changes due to displacements in the covers and along the edges. The 
derivation of the change in volume in terms of u, v, w, a, and |3 
is presented in this appendix. 

The derivation is carried out using the formula from calculus for 
determining the volume of a three-dimensional body by integrating over 
the surface enclosing the body. (See ref. 1. ) Thus,, the volume of the 
deformed plate Is given by 


V = - JJ £ dTi (Al) 

where r\, and £ denote rectangular Cartesian coordinates (along 

the x— y y-, and z-axes, respectively) locating points in the surface of 
the deformed plate. The double integral in equation (Al) is a surface 
integration carried out over all external surfaces of the plate in the 
deformed state. The negative sign in equation (Al) arises because the 
surface of integration is assumed to be oriented so that an inward 
normal is positive. 

The initial volume of the plate can be represented by the integral 


h dx dy 


(A 2) 


which is simply a double integration over the undeformed middle plane 
of the plate. The change in volume AV is given by the difference 
between V and the quantity (A2) as 


rr 


AV = 


h dx dy 


(A3) 
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Evaluation of Surface Integral 

Upper cover .- The contribution of the surface integral over the 
upper cover of the plate can be written by sutstituting 

I = x + u + (x, y^ 

t) = y + v + (x,y) \ (a4) 

t = | + 


into the integral in equation (Al). Over the upper cover, then, this 
integral is 



i d| dTi = 


dU,T)) 

d(x,y) 


dx dy 


(A5) 


where the defor ma tion of the upper surface is considered to be a transfor- 
mation of coordinates from the xyz system to the system as given by 
equations (A4). The double integral on the right-hand side of equation (A5) 
is performed over the undeformed middle plane of the plate, and the quan- 
tity ^ t s the Jacobian of the transformation from the Iq- plane to 

d(x,y) 

the xy- plane. The upper cover projects onto '.he xy- plane in a negatively 
oriented region; thus, the minus sign of equation (Al) is canceled out on 
the right-hand side of equation (A5). The Jacobian is 


d(x,y) 




bx 

dy 



cix 

dy 


( 1 + u +>x)( 1 + v +,y) " u +,y v +,x 


(A 6) 


On substitution of equations (A4) and (A 6) into the integral (A5), there 
results for the contribution of the surface integral over the upper cover 
to the volume of the deformed plate 


vo h vo tr* 
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MliJlI dx dy 
8(x,y) 



+ v 


+,Y 


+ u. 


V 

-X +,y 




+ w L + w,u 


+ +,x 


+ w + v + 


y + w + u +,x v +,y 


w + u +, y v +, x 


dx dy 

(AT) 


Lower cover.- For the lower cover of the plate the quantities 4, 

T), and £ are expressed in terms of displacements and the xyz coordinates 

as follows: 


4 = x + u_(x,y) "8 
11 = y + v_(x,y) V 

i - - | + v_(x,y)^ 


(A8) 


The surface integral over the lower cover becomes 



4 



8(4, n) 

8(x,y) 


dx dy 


(A9) 


where the minus sign must be retained because the lower cover projects 
onto the xy-plane in a positively oriented region. The Jacobian is now 

written 


8(4,4) 

8(x,y) 


8 | 

84 

dx 

1 

8 y 

dr| 

84 

dx 

8y 

( i + 

u . 
“ ) - 



+ V 



u - v v - x 
y y y x 


(A10) 


and the contribution of the integral over the lower cover to the volume 
of the deformed plate is 
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i 


d(g, Tj) 

d(x,y) 


dx dy 



+ v 


+ u 




u -,y v - 



- v_ - v_u_ ;X - w.v. ;y - V_u_^ y 


+ v - n -,y v -,x 


dx dy 


(All) 


• - It is assumed that the drop cords along the edge of the plate 
remain straight as they are displaced. Thus, the final coordinates of 
points along the edge of the plate are given by the following functions: 


I 

il 

i 


u.(x,y) - u_(x,y) 

X + — 

h 

+ v +(x,y) - v_(x,y) 
h 

v +(x,y) - w_(x,y) 

2 4 - 

h 


Z 


z 


z 


u + (x,y) + u_(x,y) 


A 


v+(x.y) + v_(x,y) 


w +( x y) + w_(x,y) 


J 


(A12) 


In these equations x and y are on the boundary of the middle plane 
of the plate and can be thought of as functions of s, the coordinate 
along this boundary. The surface integral over the plate edge can be 



£ d| di) = - 



i 


Miiii) 

d(z,s) 


dz ds 


(A13) 


The normal vector to the plate edge projects onto the z, s surface in a 
positive region; hence, the minus sign is retained in equation (A13) 
The Jacobian in equation (A13) is ' 


VO H VO (H 
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as as 

az as 

St] 5t| 

Sz Ss 

u + ~~ u - S|_ dx a^ dy 

h ax ds ay ds 


v + ~ v - a^ dx ) br\ dy 
h ax ds ay ds 



When the last of equations (A12) and equation (Al4) are substituted into 
equation (A 13 ), the contribution to the volume of the surface integral 
over the plate edges becomes 
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Total volume.- The total volume in the deformed plate is, of course, 
simply the sum of equations (A7), (All), and (A15). When equations (4) 
are used to express the deformed volume in terms of the overall plate 
displacement and rotation quantities, there is obtained 



V = / / < ( 1 + e)h 


1 + u, x + v >y + u^ x v >y - u ; .,v x + -(ha) >x (hp) >y 


+ w 


- j^( * ia ) , y ( h3 ) , x 

+ (ha)^ x v^ y - (ha) >y v^ x 
+ vh(av jX - pu >x ) - whp 


( ha) f x + (hp), y + (hJ), y u, x - (hp) >x u, y 

3 


> dx dy 



(1 + e)^(ap jX - pa x ) 


dx 

ds 


(1 + e )^ 2 ( a ^, y “ ^y) + wh Ky 


- Pu^y) + wha 


ds [ 


(ai6) 


Simplified Expression for Change in Volume 

Since a linear theory is sought, any terms in equation (Al6) which 
are of a degree higher than two in the plate displacements and rotations 
or their derivatives may be neglected. Higher degree terms lead to 
nonlinearities in the differential equations which result from the appli- 
cation of the variational technique. With this simplification the change 
in volume due to displacements in the inflatable plate becomes 


AV = Jj hU ;X + hv f y + hu fX v jy - hu jy v jX + i(ha) )X (hp) ( y - j(ha) >y (h3) ( x 

V 


+ he + w(ha) + w(hp) 


— (ap, y - Pa y ) + wha 


dx dy - j— (aP^ - P^x) - wh P 


dx 

ds 


*3L 

ds 


ds 


(A17) 


When the relations 


0\ 1 — 1 O ' \ 
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dx _ dy 
ds ~ dN 


dy dx 

ds " SnJ 


(A18) 


are utilized, Gauss' theorem in the form given in equation ( 27 ) can be 
applied to the line integral in equation (A17). The final result for 
AV, after substitution for e in terms of a and 3, is 


AV 


If ' 


u > x + Y yJ + U ,X V , y - u , y v ,x + 12 ( a ,xP,y - a , yP, x) 



§£ . 


aw 




dx dy 


(A19) 
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APPENDIX B 

ESTIMATION OF VALIDITY OF LINEAR THEORY 


A simple problem is discussed in this appendix in order to provide 
some idea of the range of parameters for which linear theory is valid. 

The results of a nonlinear analysis of the problem are compared with the 
results of the linear analysis. The problem considered is a long 
plate of constant depth which, after inflation, is simply supported on 
its long edges in such a manner that these edges cannot move. For 
simplicity the plate is assumed to be loaded by a uniform lateral load 
equally divided between the upper and lower covers. For the linear solu- 
tion, the manner in which the loading is divided between the covers is 
immaterial. For a nonlinear solution, on the other hand, if the loading 
is not equally divided between the covers, an additional nonlinear term 
arises in the theory. A cross section of the plate under consideration 
is shown in figure 3. The x-axis is assumed bo be normal to the long, 
simply supported edges, and the distance between these edges is L. 

For this problem the potential energy in varied form can be written 
as follows when nonlinear terms are included: 


511 = 


0 L 


N. 


:S( U ,X + \ + y e %) + + e ,x w ,x) 


- P h6 ( U ,X + e - aw ,x + eu ,x) - < ls 8w 


dx 


(Bl) 


where 


N v = ph + N, 


N v = 2An u 


M v 




+ i ~ e 2 y 


11 u ,x + 2 


- A ll( a jX + e ,x w ,x) 


e = - 
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(B2) 
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When the variation of equation (Bl) is performed and appropriate integra- 
tions hy parts are carried out, there results the following system of 
nonlinear differential equations: 


N x,x + P ha ^,x = 0 

N x w ,xx + N x, x w ,x + Ph(a + - (Va x )^ + ^ = 0 

w „ h 2 2 ph3 2 /„ h 2 2\ ph5 / 2 'l 

Hc,X - N X 4 “Sx ^ “x + ( N x — a ,x a ) x + — ( a ,x a ; 

- ( M x w ,x a ) >x + M x w ,x a ,x - P h ( a + v ,x) - P^X* = 0 

and the following boundary conditions at x = 0 and L: 
u = w = 0 


(B3) 


J 


? 3 

h _ ,2 , P h 2 


Mx - MxW >x a + N x — a x a^ + — a x a5 = Oj 


(B4) 


This system of equations can be simplified (partially linearized) 
by assuming that a and its derivatives are small compared with w 
and its derivatives and by neglecting terms of second degree or higher 
in a and terms which contain products of a times w or u or their 
derivatives. For practical inflatable -plate proportions, this assumption 
appears justified in this problem. When the equations are simplified in 
this manner there results the differential equations 


^x w , xx + ph(a + w^ x ) ^ + q z 
Mx,x ' P h ( a + w ,x) = 0 

and the boundary conditions at x = 0 and L 

u = w = M x = 0 


° 

i 

j 


(B5) 


(B6) 
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where 

"x - 2 A ll(",x * \ «f, 

IA - ^ A 

“ 2 ^H a ^ x 

An exact solution to this system of equations can be obtained. 

The result is in terms of rather involved tr*Jiscendental expressions, 
however, and it is difficult to find a relationship between the nonlinear 
and linear solutions. For the purposes of this appendix a simple approxi- 
mate solution using the Galerkin method suffices. 

Equations (B5) can be written 


- Cl 

(B8a) 

2A ll C l w ,xx + ph(a + v,x) )X = 

(B8b) 

2 

Y A ll a ,xx - P h ( a + w ,x) = 0 

(b8c) 



where C ^ is a constant and the boundary conditions (B6) at x = 0 
and L are 

u=w=a v =0 (B9) 

y A 

The following assumptions for w and a satisfy the last two of equa- 
tions (B9) exactly: 


w 


. irx 
W-. sin — 
L 


(B10a) 


a 


a-^ cos 


JtX 

L 


(BlOb) 


When equations (BIO) are substituted into eqi.ation (B&b) and the 

resulting equation is multiplied by sin — and integrated over x 

L 

from 0 to L, there results 


VO P vo 
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+ f (BU) 


From equation (B8c) there is obtained. 


a l = 


jt 

" V 1 L 


hf, A ll/ jt \ g 
2 ph \L/ 


(B12) 


The constant C-^ is obtained by substituting equation (BlOa) into 

equation (B8a), integrating to find u, and using the boundary conditions 
u = 0 at x = 0 and x = L. The result is 


’ til 


(B13) 


When equations (B12) and (B13) are substituted into equation (Bll), some 
algebraic manipulations lead to 




J ^ l 2 J + l 

J * P 3 


(Bill) 


A .2 

where J = -iij— | . Calculations showed very little difference between 
this approx im ate result and the exact solution. The linear solution is 

obtained when (J + l)(^j is neglected compared with unity. It is 

easy to see, for instance, that the linear solution is within 5 percent 
of the nonlinear solution provided that 


(J + 1) 


< 0.05 


This result is strictly applicable only to the problem under considera- 
tion in this appendix but does give a rough indication of what might 
happen in other problems. 
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(a) Plate having linear taper in depth and arbitrary shape in planform. 


h 



(b) Enlarged view of plate element showing taper angles. 
Figure 1.- Inflatable plate with basic coordinate systems. 






